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The properties of deterministic LOCC transformations of 
three qubit pure states are studied. We show that the set 
of states in the GHZ class breaks into an infinite number 
of disjoint classes under this type of transformation. These 
classes are characterized by the value of a quantity that is 
invariant under these transformations, and is defined in terms 
of the coefficients of a particular canonical form in which only 
states in the GHZ class can be expressed. This invariant also 
imposes a strong constraint on any POVM that is part of 
a deterministic protocol. We also consider a transformation 
generated by a local 2-outcome POVM and study under what 
conditions it is deterministic, i.e., both outcomes belong to the 
same orbit. We prove that for real states it is always possible 
to find such a POVM and we discuss analytical and numerical 
evidence that suggests that this result also holds for complex 
states. We study the transformation generated in the space 
of orbits when one or more parties apply several deterministic 
POVMs in succession and use these results to give a complete 
characterization of the real states that can be obtained from 
the GHZ state with probability 1. 



I. INTRODUCTION 

A very important part of the study of the entangled 
states of spatially separated systems, is the study of the 
transformations that are possible when using only local 
operations and classical communication (LOCC), since it 
allows us to classify entangled states and it can be used 
as one way of quantifying this resource. Two states that 
are related by local unitary transformations are consid- 
ered equivalent as far as entanglement is concerned, since 
both states can be obtained from each other and local 
operations cannot increase entanglement. The action of 
the group of local unitaries breaks the space of states into 
orbits 0| . Then, to transform a pure state into another 
state in a different orbit by local operations, we need to 
allow each party to apply a local generalized measure- 
ment, i.e., a POVM , on her part of the state. 

For bipartite pure states, the problem of deterministi- 
cally transforming a state into another has been solved 
by Nielsen who gave necessary and sufficient condi- 
tions for a given transformation to be achievable with 
probability 1 . Later Vidal || extended this result by cal- 
culating the maximal probability of success of any LOCC 
transformation of bipartite pure states. For more than 
two parties, this problem is still unsolved. The bipartite 
case seems to be very special due to the existence of the 
Schmidt decomposition. Any pure bipartite state can be 



transformed by applying local unitaries into a state of 
the form 



(1) 



where the Xi are positive real numbers, \ii) = \i)a ® 
and {\i}} are orthonormal vectors on each subsystem. 
This greatly simplifies the analysis of LOCC transforma- 
tions: it gives a canonical expression for states in a given 
orbit, and allows the reduction of an arbitrary LOCC 
protocol to a protocol in which one party applies local 
unitaries and local POVMs, and the other party only 
has to apply a local unitary, conditional on the results 
obtained by the first party |Q. For multipartite states 
with three parties or more there is no known reduction 
of LOCC protocols. 

For a system of three qubits, several Schmidt-like de- 
compositions have been proposed |||]6), all based on the 
idea of using local unitaries to get rid of as many coeffi- 
cients as possible. One interesting property that emerges 
from these decompositions is that in general it is not pos- 
sible to make all the coefficients real. In particular there 
are states that have at least one coefficient that is com- 
plex for any local basis, and this has as a consequence 
that these states are not locally unitarily equivalent to 
their conjugates (the states obtained by taking the com- 
plex conjugate of the coefficients). This contrasts with 
the bipartite case in which, since the Schmidt decompo- 
sition has only real coefficients, every state is in the same 
orbit as its conjugate. 

A POVM applied to a state has, in general, outcomes 
that belong to different orbits. However, a protocol that 
transforms a state into another with probability I, has to 
include at least one POVM for which all outcomes are in 
the same orbit. For instance, this has to be the case for 
the last POVM of the protocol: if its outcomes are not 
in the same orbit, then the protocol has not achieved the 
transformation with probability 1. We will call a POVM 
with this property a deterministic POVM, because we 
can use such a POVM and suitable local unitaries, to 
obtain any state in the orbit of the outcomes with prob- 
ability I, attaining a deterministic transformation. Since 
any local POVM can be replaced by a sequence of 2- 
outcome POVMs, it is then interesting to study the case 
of a deterministic 2-outcome POVM. 

In this paper we will study some properties of deter- 
ministic LOCC protocols and deterministic POVMs ap- 
plied to 3-qubit pure states. We will only be interested 
in transformations between states that have genuine tri- 
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partite entanglement (i.e., all three reduced density ma- 
trices have rank 2), since other cases can be reduced to 
the bipartite case. In Section II, we prove that a certain 
function of the states is invariant under any deterministic 
LOCC protocol and show that this imposes a constraint 
on the local POVMs that can be a part of a deterministic 
transformation. We also show that the set of states in the 
GHZ class breaks into an infinite number of disjoint sub- 
classes under this particular type of transformation. In 
Section III, we study the particular case of a 2-outcome 
deterministic POVM, and discuss what are the conditions 
for its existence. We prove that such a POVM can always 
be found for real states, and present some evidence that 
the same situation holds for complex states. In Section 
IV we analyze the transformation in the space of orbits. 
In Section V we study the case of the GHZ state and give 
a complete characterization of all the states with real co- 
efficients that can be obtained deterministically from it. 
Finally, the conclusions are presented in Section VI. 

II. GENERAL PROPERTIES OF LOCC 
TRANSFORMATIONS OF 3-QUBIT STATES 

Pure states of three qubits with 3-particle entangle- 
ment are divided in two inequivalent classes: the GHZ 
class and the W class [Q. They have the property that 
any local POVM applied to a state in a given class, can 
only have as outcomes states in the same class. In par- 
ticular, states in the W class can always be transformed 
by local unitary operations, into a state with real coef- 
ficients. In this paper we will call a state "real" if it is 
locally unitarily equivalent (LUeq) to a state with real 
coefficients. States in the GHZ class can be either real or 
complex. 

Any state in the GHZ class is LUeq to a state of the 
form 0| 

\ 1 p)= f x\00Q)+^\ip A )\ip B )\cp c ), (2) 

where /i > v > are real numbers, 7 G [0, 27r) and 
\ip x ) = cos6 x \0) + sin 5x11) with 5 X € (0, |] and X = 
A, B, C. We will assume that the state is normalized, 
so only five of the six parameters in (||) are independent. 
If we write \ip) = \fj,) + \v) where Jh) and \v) correspond 
to the first and second term in (0) respectively, we can 
construct the invariant 

f2(|V>)) = (fJ-W) = fj.ve n cos 5a cos 5b cos 5c- (3) 

If /.t = i>, the sign of the phase 7 is not well defined, since 
in this case there is an ambiguity with respect to which 
product state in (|^) is and hence we can interchange 
and \u) by local unitaries, and transform the state into 
its conjugate, which changes the sign of 7. As shown in 
this means that the state is real, although we need 
to use complex coefficients if we want to write it in the 
particular form given by (Q). Aside from this ambiguity, 



this decomposition is unique. If /i > v then the state 
is complex if and only if Im(f2( | tp))) 7^ 0. If Im(il(|-0))) = 

0, then either 7 is equal to or it (and in both cases all 
the coefficients are real, so the state is real), or 5x = f 
for some X. If this is the case, we can get rid of the phase 
by applying the local unitary 

to party X, which makes all the coefficients real. 

Let {Ai},i — 1, . . . ,n represent a local POVM applied 

by Alice. If we apply it to a state l^) we can write the 

_ 1 

normalized outcomes as \(j>i) = q i 2 Ai ® 1 <S> llV 7 )) where 
qi = {tl>\A\Ai ® 1 ® l|^i) is the probability of outcome 

1. Let's consider the case in which none of the opera- 
tors Ai corresponds to a projective measurement (i.e., 
they all have rank 2). If we apply this POVM to a state 
with genuine tripartite entanglement, all the outcomes 
will still have 3-particle entanglement. To understand 
why this is true, suppose that there is an operator Aj 
of the POVM such that its corresponding outcome \<pj) 
has no 3-particle entanglement. Then \tfij) has to be the 
product of a pure state of one of the parties, let's call it 
X, and a pure state (possibly entangled) of the remaining 
two parties, so party X is completely disentangled from 
the other two. Since we are assuming that Aj is invert- 
ible (it is a rank two, 2 by 2 matrix), we can construct 

a local POVM with operators {At 1 , ^1 - {Aj 1 ^ A' 1 } 

that when applied to \4>j) has at least one outcome that 
has 3-particle entanglement (the one corresponding to 
Aj 1 (gi 1 (g) l\(f>j}), that occurs with nonzero probability 

(because (Aj 1 )^ 1 Aj 1 also has rank two). Then we would 
have a protocol that with finite probability and only ap- 
plying local operations, allows us to create entanglement 
between party X and the other two, starting from a state 
in which party X was disentangled, and this is clearly not 
possible. 

Let's consider a state in the GHZ class and let Alice 
apply a local POVM to it. Then all the outcomes \<fii) 
have to be in the GHZ class too, so we know that we can 
apply local unitaries to them such that 

(U A (i) ® U B(i} <g> Uc{i))\^i) = \im) + Wi), (5) 

where 

=Mi|000) 

\vi) =v i e* li \y A ^))\(p B {i))\vc(i))- (6) 

Since \ni),\vi) and are product states, and the 

action of the POVM and any local unitaries is still local, 
for every outcome i we must have either 

VOilVi) = (U A (i) ® U B(i) U c{ i))(Ai ®1® l)\n) 
VQiWi) = ( u A(i) <8> U B {i) <8> Uc^Ai 1® l)\v), (7) 

or 
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= ® %W ® U C (i))(Ai 1® l)]ju>. (8) 

To decide which one is the case, we note that decom- 
position @ requires that /ij > Vi, and /^i,^ are the 
norms of the states \m) and respectively. Then, if 
(/il^A-Ai ® 1 (g) > (^|A]Ai ® 1 (g) l|z/) (which is equiv- 
alent to > we have that (]?]) must hold. 
Otherwise, (||) holds. Using Yli-^-t^-i — 1, we can then 
write 

Re(fi(|V») - <A*|v) + (v|A*> = ^%((^l^> + Mw» 

= £>Ite(fi(|&))). (9) 

This result is due to Vidal |)) . It puts a strong constraint 
on deterministic LOCC protocols, as we show in the fol- 
lowing theorem. 

Theorem 1 Let \ip) and |£) be two states in the GHZ 
class and assume there is a LOCC protocol that trans- 
forms \ip) into |£) with probability 1. Then, 

Re(n(|V»)=Re(Q(|0)), (10) 

i.e., the quantity Re(f2) is invariant under deterministic 
LOCC transformations. Furthermore, it must be invari- 
ant for every local POVM in the protocol, that is, if the 
POVM is applied to a state \x) and has outcomes \4>i), 
then 

Re(Q(| X »)=Re(Q(|&»), (11) 

for all i. 

Proof: The most general LOCC protocol is a sequence 
of local unitaries, local POVMs and classical communi- 
cation between all the parties. Local unitaries cannot 
change Re(f2) because is an invariant of the or- 

bit. Thus, it can only be changed by applying POVMs. 
Consider the first POVM of the protocol, that takes 
the state \ip) into one of its possible outcomes \<fri), 
each occurring with probability qi. Then, according 
to equation (||) (and because g, > 0), either all out- 
comes \4>i) satisfy Re(Q(|0j))) = Re(Sl(|?/>))) °r there 
are at least two outcomes \4>\) and \4>2) that satisfy 
Re(fi(|0i))) < Re(fi(|V>))) < Re(fi(|0 2 ))). :t is eas y to 
see that in the latter case, at any stage in the protocol, 
we will have two outcomes \<pj) and \<pk) that will satisfy 
Re(fi(|<^-))) < Re(Q(|^fe))). This will be true in partic- 
ular for the last stage of the protocol. But that would 
mean that \cj)j) and \<f>k) are in different orbits (because 
is invariant under local unitaries), and that contra- 
dicts the fact that the protocol is deterministic. Thus, 
the only possibility is that all the outcomes of the first 
POVM have the same value of Re(S7). We can apply ex- 
actly the same reasoning to all the POVMs in the proto- 
col, and then conclude that all the final outcomes satisfy 



Re(f2(|0i))) = Rc(fl(\ip))) . Since this is a deterministic 
protocol that transforms \ip) into then all these out- 
comes should be in the same orbit as and so we have 

Re(O(|0)) = Re(n(|&») = Re(fi(W)).n 

This theorem tells us that under deterministic LOCC 
transformations the class of GHZ states breaks into an 
infinite number of subclasses that are labeled by the real 
part of the complex invariant fl. Two states in different 
subclasses cannot be transformed one into the other with 
probability 1 by means of local operations and classical 
communication. From equation (|3|) and from the range 
of the parameters, we see that the set of these subclasses 
is isomorphic to the open segment (— ^, | ). The subclass 
that contains the GHZ state \GHZ) = ^(|000) + 1111>), 
corresponds to the center of the segment, and it is defined 
by Re(Sl) = 0. Note that all subclasses contain both real 
and complex states. 

This result gives a broad description of how a state can 
be transformed in the space of orbits with probability 
1. Tighter constraints can be obtained from studying 
the behavior of the entanglement monotones |]lo|] , which 
usually introduce some necessary conditions that must 
be satisfied in order for a transformation to be possible 
to be implemented locally. To find sufficient conditions 
we have to be able to show that a protocol exists that 
accomplishes the transformation. A first step in that 
direction is to study deterministic POVMs. 

III. DETERMINISTIC 2-OUTCOME POVM 

In this section we will study under what conditions a 
2-outcome POVM is a deterministic POVM (i.e., both 
outcomes are in the same orbit). A general 3-qubit state 
can be written 

l 

w>= E M^fc)- (is) 

ijk=0 

Following , we can define matrices To and T\ , where 
(Ti)jk = tijk- (13) 

The group of Local Unitary (LU) transformations of 
three qubits is locally isomorphic (i.e., has the same Lie 
algebra) to U(l) x [SU(2)] 3 . Under a LU transformation 
performed only by Bob and Charlie with matrices Ub 
and Uc, the matrices T transform according to 

T, -» UbWc, (14) 

while if the transformation is performed by Alice, we have 

T\ — > uf T + UyyTx, (15) 

where uf^ are the matrix elements of U a- 
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We know |T2|,|l3[ that the orbits of 3-qubit states can be 
parametrized with 5 continuous invariants plus a discrete 
invariant, since in general a 3-qubit state is not LUeq to 
its complex conjugate. There are many ways of choosing 
these invariants In this paper we will use the 

following set 



h 
h 
h 
h 
h 



tkijt mi jt rnpq t kpq — trp A 

ijkmpq 

y ] tikjtimjtpmqtpkq — trp B 
ijkmpq 

y ] tijkXijmtpqmXpqk = ^ r Pc 
ijkmpq 

^ ^ ^ijk^lmn^opq^rst^il^or^-jm^-ps^kq^nt\ 
ijklmnopqrst 

y , tijktii m tnlotpj tpqmt n qki (16) 
ijklmnopq 



where tij is the antisymmetric symbol and all the indices 
are summed from to 1. I4 is the 3-tangle introduced in 
JT]| . As shown in |L2| these 5 invariants are algebraically 
independent. However, since they are all real and invari- 
ant under complex conjugation of the coefficients Ujk, 
they cannot distinguish between a state and its conju- 
gate. To fix this ambiguity we use the complex invariant 



I& — y ' tiijik 1 ti 2 j2k2^i3j3k3^i4j4k4ti 5 j 5 k 5 ti 6 j 6 k e X 
iljlki 

iijik-3 i-ijiki l i3j4k s z i 4 j3ki injafo iajs^e ' V ' 

where again all indices are summed from to I. To com- 
pletely specify an orbit we need the value of I\ through 
I5 plus the sign of the imaginary part of I§. It is worth 
noting that 1 — 1\ , , 1 — 13 and I4 are decreasing en- 

tanglement monotones, while I§ is not an entanglement 
monotone fll5[ . 

We will consider the case of a 2-outcome POVM ap- 
plied by Alice on a pure state \ip) of three qubits. The 
most general POVM is given by the operators A and 
A\, where 



- l " = 1 " 1 f Vv ]U 



1 r '! J-y ]1 " 



(18) 



where Vq, Vi and U are unitary matrices JT6| , and 
< £, y<l. It is easy to see that they satisfy 
AqAq + A\Ai — 1, where 1 is the identity matrix. When 
we apply this POVM to a state \ip), we obtain two out- 
comes \4>q) and \4>i) given by 



\4>i) = —={Ai <g> 1 <g> 1|V)) » = 0,1, (19) 

V 1i 



where qi is the probability of outcome i. From ( |18[ ) and 
( |l9| ) we can see that the action of this POVM on \ip) 
is equivalent to applying a unitary tranformation first 
given by U, applying a diagonal and real POVM and 
finally applying a unitary Vi conditional on the outcome 
of the POVM. This last local unitary cannot change the 
orbit of the outcome \4>i)- Since we are considering two 
states in the same orbit to be equivalent, we can take this 
unitary to be the identity without loss of generality. 

Let us consider first the case in which U = 1. Then 
both elements of the POVM reduce to real and diagonal 
matrices 



£0 



yfx 

Vv 



E 1 



VT^x 



VT^v 



(20) 



From now on, we will take < x, y < 1, since when 
x or y are equal to zero or one, the POVM becomes 
a projective measurement, which destroys three particle 
entanglement. We can write explicit expressions for both 
outcomes of the POVM 

I </>()) = -4= ^(Vx kjk\0jk) + y/y t Xjk \ljk)) 

|0i> - ^^(%/T^ tojfclOjfc) + V^~V hjk\ljk)). (21) 
V 91 jk 

Now we calculate the invariants Ii through I 5 for |<^ ) 
as a function of x and y 



h(x,y) = 



h(x,y) = 



2 a 2 + 2xy Tt[T tI] Tr[TiT^] + y 2 b 2 

(ax + by) 2 
2 F Q + 2xy Tr[T TlT x Tl] + y 2 F 1 



(ax + by) 2 



h(x,y) = 



2 F a + 2xy TrlToTfaT*] + y 2 F l 



(ax + by) 2 



h(x,y) 



x y h(\ip)) 

(ax + by) 2 



T , v x 3 G 00 + 3x 2 y Got + 3xy 2 G w + y 3 G n 

h{x,y) = -, —TT5 , (22) 

(ax + by) 6 

where the matrices Ti are as defined in jl3]), 
a = Tr[T T$], b = Tr[T\T^], a + b = 1 for a normal- 
ized |t/>), F t = Tr[(TiT}) 2 ] and Gy = Tr [TiTj Ti T} T 3 t} ] . 
The invariants for \(j>i) are obtained from (|22j) by replac- 
ing x by 1 — x and y by 1 — y. For the two outcomes to 
be in the same orbit, we need the five invariants to take 
the same values for both states, i.e., 



h(x,y) =Ii(l-x, l-y) i=l,...,5. 



(23) 



If these conditions are satisfied, then either \<fio) is LUeq 
to \(j>i), or \(po) is LUeq to \4>i)* ■ To determine which one 
is the case, we need to calculate the sign of the imagi- 
nary part of the complex invariant Iq. For now, let us 
concentrate on the equations in (E3|). These equations 
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have a common solution with < x, y < 1 if and only if 
the following conditions are satisfied (see appendix) 



= cos a Tq + sin a T\ 



a 2 Tr\(T x Tlf\ = b 2 Tr[(T Q T^) 2 } 

a Tr[Ti Tj 7\ t\ T t\ ] = b TrpbT^ToltfTiT,}] 
a 2 x{l-x) = b 2 y{l~y). 



(24) 
(25) 
(26) 



Furthermore, the solution satisfies is ( | ) ) < IsQtp)). 
This is worth noting because I5 is not an entanglement 
monotone, but behaves monotonically under this partic- 
ular class of POVMs. Equations (fj) and (f|) are real 
valued polynomial constraints on the coefficients of the 
state, and in general are not satisfied for an arbitrary 
state. From (|lj) and (|l5|) we can see that these con- 
straints are invariant under LU transformations applied 
by Bob and Charlie, while they are not invariant under 
local unitaries by Alice. Equation (|2^) is a constraint on 
the parameters of the POVM that depends on the state 
we are transforming. 

Now let U be any unitary matrix, so our POVM takes 
the form {E U, E X U}, with E ,E X given by (j20j) . This 
is equivalent to applying the local unitary U to Alice's 
part of the state, followed by a diagonal POVM, and 
we know the conditions that need to be satisfied in this 
last stage. So we can reduce the problem to finding a 
local unitary performed by Alice that would transform 
the original state into a state that satisfies (|3|) and 
). Then we can choose a POVM that satisfies (26), 



(| 

where now a and b are calculated using the coefficients 
of the transformed state U <£> 1 ® 1|V')- We will consider 
the cases of real and complex states separately. 



sin a Tq + cos a T\ . 



(28) 



If we plug this into (p4|), take out a common tactor cos" a, 
introduce the variable z = tana and move all terms to 



A. Real States 



one side, we can write (24) as polynomial pi(z) of degree 
8 with real coefficients, of the form 

Pl (z) = A(l - z s ) + B(z + z 7 ) + C(z 2 - z 6 ) + 

+D{z 3 + z 5 ) = 0. (29) 

If zq is a real root of pi, then U(ao), with ao = arctan(zo) 
is the unitary matrix we are looking for. Now it's easy to 
check that pi(l) = — pi(— 1), so p x has at least one real 
root in [—1,1], which means that we can always find a 
unitary U, such that \ip') = U (g) 1 ® l\ip) satisfies (||). 
Now we can apply to a diagonal POVM that satisfies 
(ptf), and we are certain that both outcomes have the 
same values of the four invariants = 1, . . . , 4. But in 
the case of real states this is enough to completely specify 
the orbit, because since \ip) is real, so is \ip') because U 
was chosen to be real, and the outcomes of the POVM, 
\4>o) and \4>i), are also real, because the POVM itself 
is real. In this case we don't have to worry about the 
value of the complex invariant. Finally, since \<fio) and 
\4>i) are in the same orbit, we can apply local unitaries 
to transform them into any state in the same orbit. So 
the net result of this protocol is to transform any state 
in the orbit of \ip) into any state in the orbit of |<^o), with 
probability 1. The results presented so far show that for 
any real state, there is some set of orbits that can be 
reached deterministically from that state, although we 
haven't yet characterized this set. We will discuss this 
problem in Section IV. 



To characterize the orbit of a real state \ip) we only 
need four parameters instead of the five needed for an 
arbitrary state. First, note that, by our definition, any 
real state can be transformed by means of local unitary 
transformations, into a state with only real coefficients. 
Of the (at most) eight coefficients of this state, only seven 
are independent if we are considering a normalized state, 
and we can get rid of three more by applying local real 
unitary (orthogonal) transformations on each of the 3 
qubits. Since Jj,i = 1, ... ,4 are algebraically indepen- 
dent, we can use this set to parametrize the orbits of 
real states. This greatly simplifies our analysis because, 
as seen in the appendix, ( pi] ) is enough to assure that 
Ii,i = 1, ... ,4 have the same values for both outcomes 
of our POVM. So, given a real state, we need to find a 
U such that \i/>') = U ® 1 ® satisfies @. Let 

t cos a sina\ 
v ' V — sin a cos a 1 y ' 

In terms of the matrices Tj, this transformation can be 
written 



B. Complex states 

The analysis of the complex states turns out to be more 
complicated, because now we need to find U such that 
satisfies both (|i]) and (|25|). We can write any unitary 
U as 

i<t> ( eW \ f C0Sa Sina \ ( e< ^ (QM 

\ e-* 3 J \ - sina cosa ^0 e~* ) ' 1 ' 

The phase and the matrix on the left will commute with 
the diagonal matrices of the POVM, so their action is 
equivalent to applying a local unitary to the outcomes of 
the POVM. But we know that this action will not change 
the orbit of the outcome state, so we can fix them to be 
the identity. So U will take the state with matrices 
To and T\ to a state \tp') with matrices 

Tp 1 = cosa e iC T + sina e~ iC T x 

T[ = - sin a e lC T + cos a e' iC T x . (31) 

We can substitute (|3l]) into the homogeneous form of 
( g4j ) and (|25|) , again divide by cos 8 a and introduce the 
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variable z — tana, so both conditions are expressed as 
polynomials in z equal to zero, with real coefficients, of 
the form 

Pi (z) = Ml - z 8 ) + B t {z + z 7 ) + a(z 2 - z 6 ) + 

+D t (z 3 + z 5 )=0 i = l,2, (32) 

with the coefficients given by 

Ai = aoi 

Bi = bu cos(2C) + b 2l sin(2C) 

Ct = c 0l + c u cos(4C) + c 2i sin(4£) 

Di = du cos(2^) + d 2i sin(2() + d 3i cos(6C) + 

+da sin(6C), (33) 



where a i, 6^, c jiy dji 



are real valued polynomials on the 
coefficients of whose exact form can be computed 
from regrouping the terms obtained after substituting 
(|l|) into dJ) and (H). 

Finding a local unitary performed by Alice on \ tp) that 
would yield a state that satisfies (^J|) and (|||), is equiv- 
alent to finding values z and C (which parametrize the 
unitary) such that both polynomials p\ and p 2 vanish. 
We can think of ( as a parameter for these polynomials, 
and what we are looking for is a value of £ such that p\ 
and p 2 have a common real root. 

The polynomials pi have certain useful symmetries. 
First of all, because their coefficients are real, complex 
roots appear in conjugate pairs. Also, because of the 
particular symmetry of the coefficients (i.e., the coeffi- 
cient of z s is equal to minus the independent term, the 
coefficient of z 7 is equal to the coefficient of z, and so on), 
if zo is a root of pi, so is — j- (this corresponds to ct+ | 
being also a solution if so is a, and this is related to the 
fact that ( p4| ) and ( |25| ) are symmetric under the inter- 
change of and 1). Since Pi{l) = — Pi(— 1), Pi has a real 
root in [—1,1]. To simplify the problem, we can extract 
a factor z 2 + 1 from pi , so we reduce the problem to two 
polynomials of degree six, that have the same symmetry 
properties discussed above. If we introduce the variable 
w = (- — z) we can further reduce the two polynomials 
of degree six to two polynomials gi(w), i = 1, 2 of degree 
three, given by 

gi(w) = A,w 3 + B lW 2 + (C t + 2Ai)w + (A + Bi), (34) 

with the property that if w is a root of gi, the correspond- 
ing z's given by w — (- — z) (which are real if and only if 
w is real) are roots of pi. So we reduced the problem to 
finding a common real root of g\ and g 2 . The resultant 
Ji7t of the two polynomials g\ and g 2 is a function of C, 
and takes the form 



Ree(ffi, saXC) = I> fe cos t( 2 + 4fc Kl + Sk sin t( 2 + 4fc KD> 

(35) 



where Tk and Sk are polynomials on the coefficients of ■ 
We can see that this resultant vanishes several times in 
[0, 27r], which is the range of (, and this is useful because 
the resultant of two polynomials vanishes if and only if 
they have a common factor. This falls short of saying 
that we can find £ such that <?i and g 2 have a common 
real root, because there is in principle the possibility that 
the common factor is a polynomial of degree 2 irreducible 
over the real numbers, so g\ and g 2 have a common root 
but it is complex. However, after checking this procedure 
with many randomly generated states, we found that the 
common factor always corresponds to a real root. 

Let's assume that in fact, we can always find a value 
Co such that p\ and p 2 have a common real root zq. 
Then we know that if we apply U(ao, Co) <S> 1 <8> 1 (where 
Q'o = arctan(zp)) to we obtain a state |?/>') that sat- 
isfies (^J) and (p5|). Then, we can choose a POVM that 
satisfies (p6[), and we can be sure that both outcomes 
of this POVM, when applied to will have the same 
values of Ij, i = 1, . . . , 5. However, as we pointed out be- 
fore, this is still not enough to say that both outcomes 
are in the same orbit. There's still the possibility that 
they are in orbits that are conjugate to each other, since 
we are dealing with complex states, which are not LUeq 
to their conjugates. To decide which one is the case, 
we can calculate the sign of the imaginary part of Iq for 
both outcomes. Unfortunately, the expression of Ig for 
both outcomes is too complicated and it's not possible 
to extract the sign of the imaginary part analytically for 
an arbitrary state, although it is very easy to compute it 
numerically for a given state. We analysed randomly gen- 
erated states, and found that we can always find a value 
of C for which both outcomes are indeed in the same or- 
bit (although there are other values of ( for which the 
outcomes are in conjugate orbits) . We will refer to states 
with this property as gate states, since we can use them 
as a gate to leave one orbit and move to another with 
probability I. 



IV. THE TRANSFORMATION IN THE SPACE 
OF ORBITS 

We can now use the results of the previous section to 
give a characterization of the states that can be obtained 
from l^) by applying a 2-outcome deterministic POVM. 
Let us assume that the state is a gate state. We will 
also assume that a < b (if it's not, we apply a bit flip on 
Alice's qubit, which interchanges the matrices To and T\, 
and hence a and b) . We can use the invariants evaluated 
for |</>o) (given by (H||)) to characterize the orbit of the 
outcomes. These equations are homogeneous of degree 
zero in x and y, so we can write them in terms of only 
one parameter A = - 



fe=0 



fa 



X 



(a + 6A) 2 



i = l,..., A 



G 



A(/3 5 + 75 A) 
(a + bX) 3 ' 



(36) 



where 



, Tr[(T Ttf] 

Oil = 1 , 0L2 = «3 = 5 , a 4 = U 

gt 

rr[(T T,J) 3 ] 



2(Tr[T T 1 t ]Tr[TiT t ] - ab) 

Tr[{T Tlf] 



Or, 

ft = 2(Tr[T Q T^T 1 T}} - b 

a 

ft = 2(Tr[r T 1 t T 1 r o t ] - b Tr ^ T ° T ^ 

a 

ft = 3(Tr[r r 1 t Tor t T 1 T t ] - brr[(T ° T( 



t\3l 



75 = l{Tr[T x TlT x T\T,T\\ - ^K^o) 3 ] 



')■ (37) 



The range of A is [1, +oo) (when a <b), with A = 1 cor- 
responding to no transformation (Eq oc 1), so we have 
ij(A = 1) = and A = +oo corresponds to a pro- 

jective measurement (y — l,x — 0), that destroys any 
3-particle entanglement. From ( |36| ) we can see that the 
set of orbits we can reach from can be described as a 
one parameter family {Ii(X)} that corresponds to a curve 
in the space of orbits, that starts at state \ip) and ends 
on a state that has no tripartite entanglement. 

It is possible for some orbits to have more than one gate 
state. The values of the coefficients (37) will be in general 
different for different gate states. Since these coefficients 
determine the curve {Ii(X)}, we will be able to transform 
to different sets of orbits depending on which gate state 
we use. We can also reach a different family of orbits if we 
let Bob or Charlie apply a deterministic POVM instead 
of Alice. This is because the matrices are different 
for different parties, and so will give in general different 
gate states. 

If we fix the sign of the imaginary part of I§, we can 
use the invariants {Ii, i = 1, . . . , 5} as coordinates for the 
orbits. All the previous results can be summarized in 
the following picture. Every point in this space (which 
represents the orbit of some state is the starting 

point of a finite number of curves, each representing a set 
of orbits that can be obtained from with probability 
1 with a local 2-outcome POVM. 

More orbits can be reached if several rounds of deter- 
ministic POVMs are allowed. The general protocol will 
be something like this: (i) starting with the state l^), 
Alice applies a local unitary to transform it into a gate 
state; (ii) she applies a POVM on her part of the system, 
that satisfies (p6|); (iii) according to the outcome she ob- 
tains, she communicates to Bob and Charlie the state 
they are sharing after the measurement; (iv) they de- 
cide which one will apply the next POVM and repeat the 



same steps, now starting with the state A simplified 
pictorial representation of this transformation is given in 
Figure 1. 



h 




h 

FIG. 1. Transformation of states in the space of orbits. 
Full lines represent a particular transformation that takes \ip) 
to \ip') and then to Dotted lines correspond to other 

possible choices for deterministically transforming the states. 

The transformation occurs in the 5-dimensional space 
defined by the invariants ii, but for simplicity, we repre- 
sent only two of them (I4 and 1$). We start with a gate 
state |^) and we apply a deterministic 2-outcome POVM 
(with some parameter Ao), that transforms it into state 
\ip'). The line connecting and \ip'} represents all the 
orbits that can be reached from \tp) by applying a POVM 
with parameter A between 1 and Ao- The dotted lines 
originating at \tp) represent the set of orbits that can be 
reached from the same orbit, but using a gate state dif- 
ferent from \ip) (which is in the same orbit as \ip) so it's 
represented by the same point in the plot). In the actual 
space of orbits, these curves extend until they reach an 
orbit that represents a state with no 3-particle entangle- 
ment, that corresponds to the point where the POVM 
becomes a projective measurement (i.e., A = +00). For 
clarity, we are only plotting the beginning of these curves. 
After deterministically transforming l^) into |t//), the 
parties can choose again from several gate states to ap- 
ply the next POVM. This will determine which party will 
apply this POVM, because in general, a state is a gate 
state only for a particular party. In the figure, the full 
line represents a POVM that transforms into IV 1 ") ; 
while again, the dotted lines correspond to other possible 
deterministic transformations that can be applied to 
By applying many deterministic POVMs with different 
parameters, we can reach many different orbits. 
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V. TRANSFORMATION OF THE STATE 

|GHZ) = ^(|000) + | HI)) 

As an example of the use of 2-outcome determinis- 
tic POVMs, we will now study the particular case of 
the \GHZ) state. As it was mentioned in Section II, 
this state belongs to the subclass of states that satisfy 
Re(51) = 0. We will show that it can be transformed 
with probability 1, into any real state in that subclass. 

First, we need to identify the real states that satisfy 
Re(fi) = 0. Then clearly we must have that either n(\ip)) 
is zero or pure imaginary. In the former case, this means 
that (pb\v) = 0, and then decomposition (|J) takes one of 
the following forms: 



M |ooo) + Hi)b)|^') 

M |000) + %>)|1>|^'> 
M |000) + %>)|^>|1) 



(38) 



If n(|^)) is pure imaginary, then the only case in which 
\ip) is actually a real state is the case in which /i = v 1 as 
discussed in Section II. In this case, the state takes the 
form 



i=(|000)±i|<^>"», 



(39) 



where none of the states in the second term can be equal 
to |0) or |1) (otherwise it could be transformed into a 
real state by a local unitary), and we obtain \i = 
by imposing normalization of the state. The two states 
in (39) (corresponding to the two possible signs of the 
second term) are LUeq to each other. 

Since the GHZ state is symmetric under a permutation 
of the parties, it is clear that if we find a protocol that 
transforms it into the first state in (38), then we can also 
transform it into the other two. In this section we will use 
the results of Section III to explicitly construct protocols 
that transform the GHZ state into the state 



or the state 



V2 



n\0M) + u\l)\<p)\<p') 



(|ooo>+V»l^», 



(40) 



(41) 



for all allowed values of fi,i/,\tp),\(p') and \(f"). These 
protocols will be divided into three steps. First, Char- 
lie applies a local deterministic POVM that transforms 
\GHZ) into -^_(|000) + \11)W))- Then, Bob applies an- 
other local POVM that takes the state to -^(|000> + 
\l)\ip)\<p')). Finally, Alice applies the last POVM, which 
she can choose to take the state to ^|000) + i/|l)|<^)|y>') 
or^(|000>+^")|^')). 



by 



Step 1. The Ti matrices for the GHZ state are given 



75° 




Tl = I A. 




1 

V2 



(42) 



and they have the same form for all parties. If Charlie 
applies a local unitary U on its qubit, where 



U : 



V2 ( 1 1 
2 l-ll 



(43) 



the T[ matrices for the state \tp') = 1 <£> 1 ® U\GHZ) are 



rpl 

1 n 



1/10 
1 



1/-10 
1 



(44) 



It is very easy to see that these matrices satisfy equation 
(p4|), so the state \tp') is a gate state. Thus, Charlie can 
apply a deterministic POVM to it. In particular, this 
state satisfies b' = a' = Tr-p^To] = h, so according to 
equation ( p6|) we have y = 1 — x, so Charlie can apply a 
deterministic POVM of the form 



En = 







y/T= 



1 l^jU W 



where x £ [i, 1). The normalized state corresponding to 
the outcome zero is 

|^o> = V2(l ® 1 ® Sb)|V0 

= >/2(l® 1®E U) | GHZ) 

= |00) (£W|0)) + |11) (E U\l)) 

= (0|(7 t ^ {/|0)^|00)|0'> + 

+{11^4^11)^11)11'), (46) 



where 
|0') = 



E o U\0) 



(o\WEiE u\oy- 



II') 



E U\1) 



{l\WElE U\l)i 



(47) 



are normalized states. A straightforward calculation 
shows that (0\W E f E U\0) = E f E U\l) = |, so 
we can write 



1 



|0o) = ^(|OO)|O') + |H)|l'))- 



(48) 



This state can be taken to the canonical form (g) by 
letting Charlie apply a local (real) unitary on his qubit, 
that takes the state |0') into |0), and |1') into \<p') — 
cos(5'|0) +sin5'|l). Thus, (0|y/) = (0'|1') and then we 
have 



cos 5' = (0'|1') 



2(0\U^ElE U\l) 
2a; -1. 



(49) 



We can see that for any 5' G (0, we can find x € [i, 1) 
that satisfies this equation. This means that we can 
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transform \GHZ) into -^(|Q0Q) + |H)|v'» with Prob- 
ability 1, for any \<p'). 

Step 2. In this step Bob applies a deterministic POVM 
to transform the state \<f>) = ^=(|000> + |H)|<p')) into 
^(|000) + |1)|^)|<^')). The Ti matrices for \<j)) from Bob's 
point of view, are given by 



T = 





°) 




( V2 




V o 


o ) 





75 cos5 ' 7f sin ' 5 ' 



(50) 



First, Bob applies the local unitary U from ( ff3| ) to his 
qubit, obtaining the state \4>') = 1 (£> U <£> l|</>), character- 
ized by matrices T/ given by 



Consider now the case in which Alice wants to obtain 
-^(|000) + i\(p")\(p)\(f')) from |%). In this case we can 
construct the appropriate POVM {Aq, Ai} by inspection. 
If \<p") = cos 8" 1 0) + smtf"|l), we define 



A 



J_ / 1 i cos 5'' 
~ ^2 V isin(5 " 



-4i 



J_f I -icosS" 
~ y/2 V -* sin ^" 



(56) 



It is easy to verify that they satisfy AjAo + = 1, 

and that the probabilities of both outcomes are equal to 
The normalized state that corresponds to outcome 



2 

zero is 



T'-i 



1 

cos 8' sin 8' 



r = l( -i o 

' 1 2 \ cos <5' sin 5' 



(51) 

Again, it is easy to show that T[ satisfy (p4|), so is a 
gate state. We also have that a! = Tr\T^T' n ] = i, so Bob 
can apply the POVM of equation ([45]) to his qubit and 
obtain two outcomes in the same orbit. We can apply 
the same analysis we did in Step 1 to the outcome |xo) 
of Bob's POVM, and show that 



1 



| Xo ) = _(|0)|0')|0) + |l) 



(52) 



where the normalized states |0') and |1') are also given 
by (|47j). It should be clear from Step 1 that, again, we 
can choose x and a suitable local unitary on Bob's qubit 
to transform this state into 



Ix) 



1 



V2 



(|000) + \1)\<pW)) 



(53) 



for any \tp) = cosc5|0) + sin<5|l), with 8 S (0, §]. 

Step 3. Now Alice has to choose between two local 
POVMs depending on whether she wants to obtain j40| ) 
or ([jl]). Consider first the case in which she wants to 
transform \x) into /i|000) +v\l<p(p'). The T, matrices for 
|x) from Alice's point of view are 



J_ / 1 

V2\0 



_ __ 1 / cos 8 cos 8' cos 8 sin 8' 



(54) 



These matrices already satisfy equation ( p4| ) and since 
a = Tt[TqTq] = |, Alice can apply the deterministic 
POVM given by (f45|). The state corresponding to out- 
come zero is 



10 = V^looo) + VT - x\i<p<p') 

= /Lt|000) + z/|l W ')> 



(55) 



where we set /x = y/x and v 
we have fj, > v. The state in 



1 — x. Since x S [ jr, 1 ), 
I) is the same as in W~ 



V2 



(|000)+i|^')b)l^)) 5 



(57) 



while the one corresponding to outcome 1 is just the com- 
plex conjugate of (|57|). But we know that these two states 
are actually in the same orbit, so we can transform out- 
come 1 into ( |57| ) by local unitaries, so we obtain (57) with 
probability 1. The state in (|57|) is the same as in pi]). 
This concludes the protocol. 

Note that all three steps involve only local unitaries 
and deterministic POVMs, so these protocols allow Al- 
ice, Bob and Charlie to transform the GHZ state into any 
other real state that belongs to the subclass defined by 
Re(£7) = with probability 1, using only local operations 
and classical communication. This is then a complete 
characterization of the real states that can be obtained 
from the GHZ state, since by Theorem 1 we know that 
we cannot reach real states that belong to a different sub- 
class. It is interesting to note that it does not seem to be 
that easy to find a deterministic protocol to transform 
the GHZ state into any complex state in the same sub- 
class. Whether this is actually possible is still an open 
question. 



VI. SUMMARY AND CONCLUSIONS 

In this paper, we studied the properties of determin- 
istic LOCC transformations of 3-qubit pure states with 
tripartite entanglement. We showed that the set of states 
in the GHZ class breaks into an infinite number of dis- 
joint subclasses, characterized by the real part of a com- 
plex function £l(\ip)). Two states that belong to different 
subclasses cannot be transformed one into the other with 
probability one, by means of local operations and classi- 
cal communication. This quantity is not only invariant 
under deterministic transformations, but it also must be 
conserved by any local POVM that is part of a deter- 
ministic protocol. This imposes a strong constraint on 
the POVMs that can be used for deterministically trans- 
forming a given state. 

It is interesting to point out that the invariance of 
Re(f2) under deterministic LOCC transformations (and 
its invariance under any local POVM that is part of such 
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a transformation) , follows from the invariance of under 
local unitaries and the very particular form of equation 
(0) . In the language of entanglement monotones, we can 
say that Re(f2) is both an increasing and decreasing en- 
tanglement monotone. Any function of the states that 
is invariant under local unitaries and satisfies an equa- 
tion like (JsJ) for an arbitrary local POVM, will be in- 
variant under deterministic LOCC protocols, and hence 
will break the set of states into inequivalent classes that 
will be labeled by that function. This will be true even in 
the multipartite case, so identifying quantities with these 
properties could be very useful in the study of determin- 
istic transformations of entanglement. 

We also discussed the case of a deterministic 2-outcome 
POVM. We showed that for this POVM to exist, both the 
state and the parameters of the POVM have to satisfy 
certain polynomial conditions. In particular the coeffi- 
cients of the state have to satisfy two polynomial con- 
straints. To be able to apply a deterministic POVM to 
a given state, we need to find a local unitary that will 
transform our original state into another state that sat- 
isfies the two constraints. For real states, the problem 
actually simplifies and only one constraint has to be sat- 
isfied. In this case, it was proven in general that the nec- 
essary local unitary could be found, allowing us to apply 
a local 2-outcome POVM that would send the state to 
some other orbit with probability 1. For complex states 
we found some analytical evidence that the unitary could 
be found, but a rigorous proof of this fact is still an open 
problem. However, it is important to stress that of all 
random numerical examples analysed, the algorithm dis- 
cussed in Section II never failed to find a gate state for 
complex states. We also discussed how several rounds 
of POVMs and local unitaries applied in sequence by all 
the parties allow us to reach a bigger set of orbits than 
the one we get from only one POVM. There is a lot of 
freedom in choosing the order in which the parties apply 
a POVM and which POVM they choose. Although it is 
in general difficult to study this procedure analytically, 
in order to characterize the set of states that can be ob- 
tained from \ip) (except for states with high symmetry 
like the GHZ state), a numerical analysis is easy to im- 
plement, and can be used to study general properties of 
this set, that could help us to have a better understand- 
ing of deterministic transformations. 

Finally, we combined the two main results of this pa- 
per to give a complete characterization of the real states 
that can be obtained from the GHZ state with probability 
1. First we used the results of Section II to character- 
ize the subclass of states that could in principle be ob- 
tained deterministically from it, and then we constructed 
an explicit protocol that allows the three parties to trans- 
form the GHZ state into any real state in that subclass. 
Finding a protocol to transform it to a complex state in 
the same subclass does not seem to be as easy, and thus 
whether this transformation is possible or not is still an 
open question. 
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APPENDIX A: SOLUTION OF 

I i (x,y)=I 4 (l-x,l-y) 

We want to know under which conditions does Q25 
have a nontrivial solution (i.e., x ^= y and x,y ^ 0,1). 
We will consider only states that have 3-particle entan- 
glement, which means that a, b ^ 0, 1. First, let us note 
that we can write Ii(x,y) as 

[ax + byY 

where (Tr[T Tl]Tr[TiT^] ~ ab) ^ if \ip) has 3-particle 
entanglement. Then I\(x,y) = — x, 1 — y) has a 
solution if and only if 



xy 



(l-x)(l-y) 



(A2) 



(ax + by) 2 (a(l-x) + b(l~y)) 2 ' 

which is the same as 

a 2 x{\- x) =b 2 y{l-y). (A3) 

This also implies that I±{x, y) = I^il — x, 1 — y). Both 
I 2 and I 3 have the form 



Tf v F x 2 + F lV 2 + 2C iX y . 
(ax + byy 



(A4) 



We can use ( |A2|) to write Ii(x, y) — Ii(l—x, l—y), i = 2, 3 
as 

F + F lZ 2 F + F lW 2 



(a + bz) 2 (a + bw) 2 



where we introduced the variables 



and w 



(A5) 
(i-v) 

(l-x)- 

From (A3) we see that these variables are not indepen- 
dent, and satisfy the condition zw — (f ) 2 . Furthermore, 
both z and w are posi tive , since x and y are between 
and 1. If we expand ( |A5| ) and use the relationship be- 
tween z and w, we have 



(F b 2 - F ia 2 ){z 



(a 2 + b 2 ) 



+ 2t)=0, 



(A6) 



and since z has to be positive (and a and b are positive), 
we have the condition 
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a 2 F = b 2 F 1 , (A7) 

which is equation (p4|). 

To study the equati on h (x, y) = /s(l — x, 1 — y) we 
can assume that both (A7) and (AS) are satisfied, since 
we are looking for a simultaneous solution o f (B3I) . Let 
ji = Is(x,y). Introducing z = - and using (|22|) we can 
write 

Goo + 3G iz + 3G w z 2 + G n z 3 = fi(a + bzf, (A8) 

where Gy — Tr[TiTjTiTjTjTj], and we can expand this 
into 

(Goo - Ma 3 ) + 3(G i - \±o?b) + 3(G 1Q - fiab 2 )z 2 + (A9) 
+(G U - t ib 3 ) = 0. (A10) 

A root of this cubic polynomial represents an operator 
of a POVM for which the value of I5 for the outcome 
of that operator is /i. We are looking for two operators 
whose outcomes have the same value of I5 , but that also 
satisfy equation ( |A3[). That is the same as finding two 
roots zq and z\ of (|A9|), that satisfy the condition 



Expanding this inequality we get 
ZqZ\Z2 - (zqZi + z z 2 + Z1Z2) + (zo + z\ + z 2 ) - 1 > 0. 

(A17) 

All the symmetric polynomials on the roots of a polyno- 
mial equation can be written in terms of the coefficients 
of that polynomial, so we can rewrite this inequality as 



(Goo ^ Ma 3 ) 3 (Gqi - fia 2 b) _ (G10 



fiab 2 ) 



1 >0. 



(G u ~ /J-b 3 ) (G n - fib 3 ) (Gn - M& 3 ) 

(A18) 

Expanding this and using a + b = 1 we get 

Goo + 3 G01 + 3 do + G u > (j,. (A19) 

But the left hand side is just the value of I5 for the state 
\tjj), while /i is the value of 1$ for the transformed state 
\(j) ) (or \4>i)). So this condition is telling us that under 
a deterministic 2-outcome POVM, I5 behaves monoton- 
ically, even though it is not an entanglement monotone 
in general. 



z o 2 i - p- 



(All) 



Let Z2 be the third root of ( A9). From elementary algebra 
we know that the product of the three roots is equal 
to minus the quotient of the independent and the cubic 
coefficients, so we can write 



Z0Z1Z2 



Gqq ~ ^ a'^-fi 
Gn - lib 3 



b 3 % 

6r 



(A12) 



Using (A7) and the Cay ley-Hamilton theorem, it can be 
shown that 



Gqq 
a 3 



so (A12) reduces to 



Z0Z1Z2 



Gn 



(A13) 



(A14) 



If we want (All]) to be satisfied we need z 2 — — f . If we 
plug this into (|A9|), we find that z 2 is actually a root if 
and only if 



b G01 = aFi 



10, 



(A15) 



which is equation (p5|). There is one more detail we need 
to check. We need zq — | and z\ = j5f to be positive 
numbers, because x and y are between and 1, and only 
one of them should be greater than 1 (which can be seen 
from their explicit form in terms of x and y). We know 
that the other root z 2 = — | is negative, so the condition 
for only one of them to be greater than 1 can be written 



(zo - l)(2i - l)(«a - 1) > 0. 



(A16) 
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